The Phase Diagram of Disordered Vortices from London Langevin Simulations 
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We study the phase diagram of vortex matter in disordered type-II superconductors. We performed 
numerical simulations in the London Langevin approximation, using a new realistic representation of 
the disorder. At low magnetic fields we find a disentangled and dislocation free Bragg-glass regime. 
Increasing the field introduces disorder-driven entanglement in a discontinuous manner, leading to a 
vortex-glass phase, which subsequently melts into the vortex liquid. The obtained phase boundaries 
are in quantitative agreement with the experimental data. 

PACS numbers: 74.60. Ge, 74.25.Dw 



Fluctuations in the High Temperature Superconduc- 
tors (HTSC) play a dominant role in vortex phenomenol- 
ogy due to the high transition temperatures Tc ~100K 
and large anisotropies e _1 ~5-500. As a result, in clean 
systems the Abrikosov lattice melts into an entangled 
vortex liquid with increasing temperature. This transi- 
tion was predicted on the basis of Lindcmann type cal- 
culations |l],|2| and is by now firmly established by Lon- 
don |[ , 3D XY j§, and Ginzburg-Landau f§ simula- 
tions. Experimentally, the latent heat found in YBCO ||] 
and the jump in the magnetization in BSSCO j7|,|| and 
YBCO H,[hJ provide thermodynamic evidence for a first 
order melting transition at low fields. 

Although disorder was first argued to eliminate the 
long range order of the Abrikosov lattice altogether [jll) , 
more recent theories propose that quasi-long-range trans- 
lational order is still preserved at low magnetic fields 
in a dislocation free "Bragg-glass" phase |I^,|r|. In- 
creasing the field effectively increases the disorder and as 
the pinning energies overcome the elastic energies, vor- 
tices entangle |T3] forming a vortex-glass [jlf|. Experi- 
mental support for this scenario is the sharp enhance- 
ment of the critical current from magnetization measure- 
ments Jl6[-^9| , the rapid destruction of the Bragg peaks in 
neutron scattering p(J , and the pronounced dips in the I- 
V curves Q . A numerical study found evidence for pairs 
of point-like disclination defects driving the Bragg-glass 
to vortex glass transition in layered systems with strong 
pins 0]. Transport measurements [§3],|4| indicate that 
the first order melting transition from the Bragg-glass to 
the vortex liquid at low fields is replaced by a continuous 
melting transition from the vortex glass to the vortex liq- 
uid at higher fields. This basic feature is reproduced by 
numerical simulations Whether the vortex glass 

and vortex liquid are separated by a phase transition or 
a crossover is presently under discussion and seems to 
depend on the screening length p7| ]. 

In this Letter we report continuum Langevin simula- 
tions for disordered vortex systems in the London ap- 
proximation, which is valid outside the 3D XY critical 
region and for fields B <0.2iJ C 2(T). This approach is 
more realistic than lattice simulations due to the absence 



of intrinsic pinning, yet it is capable of reaching satisfac- 
tory system sizes. Furthermore, by using a coarse grained 
pinning landscape, we are able to reach the weak collec- 
tive pinning regime. We measure the in-plane structure 
factor 5(Q), the vortex wandering along the field B(z), 
as well as the time overlap correlator C(t) (see below for 
the definitions of these quantities). The central result of 
our paper is the determination of the phase diagram for 
realistic YBCO parameters (Fig. 1). 




FIG. 1. The B-T phase diagram for YBCO parameters, 
consisting of a Bragg Glass, Vortex Glass, and Vortex Liq- 
uid phases. For comparison, experimental data is shown as 
dashed (Ref. Jl|) and dash-dotted (Ref. |HJ) lines. 

We start by outlining our numerical approach. The 
magnetic field is taken to be in the c-direction, perpen- 
dicular to the Cu02-layers, indexed by z. The Langevin 
equation describing the overdamped motion of vortices 
reads 



V- 



dT L^z,t) 
dt 



SH[{R v (z,t)}] 



(i) 



were [i labels the vortices with coordinates R. The damp- 
ing force on the left hand side is characterized by the 
coefficient ij that is related to the Bardeen-Stephen flux- 
flow resistivity by pes = B&q/(c 2 t]), where <&o = hc/2e 
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is the flux quantum. The Langevin force £ is correlated 
as (C2(z,t)CS(z',t')) = 2r,TS a > 3 6^S(z-z')S(t-t'), with 
a, /3=x,y. The Lorentz force from a bias current density 
J is Fl = ^oJ x z/c. The Hamiltonian H is constructed 
on the basis of London theory. Its derivative decomposes 
into three forces. The pairwise interaction force is the 
periodic extension of the in-plane London interaction, in 
Fourier space V(k) — 2ire exp(— fc 2 £ 2 )/(A~ 2 + fc 2 ), with 
line energy e = $o/(87rA) 2 , London penetration depth 
A, and coherence length £. We choose k=A/£=100. The 
bending forces are governed by the single vortex elas- 
tic constant C44 = e 2 eolog(K), with anisotropy e=l/5. 
Finally, the dimensionless pinning force on vortex fi in 
layer z for STc disorder is given by 

F£ z = -J d 2 R'V R ^ z)P (R^z) - R>(R', z) , (2) 

where the microscopic disorder potential is delta- 
correlated, (u(R, z)it(R', z')) — 7<5(R — R')<5 Zj2 ', with 
strength 7. With our conventions, the parameter 7 
equals the parameter S a in Chapter IIIC of Ref. [|). 
The convolution with the long-range vortex form factor 
p(R) = 2£ 2 / (R 2 + 2£ 2 ) generates a smooth pinning land- 
scape. 




FIG. 2. A part of a typical pinning landscape ( "pinscape" ) . 
Shown is the pinning force in the a;-direction. 

The z coordinate is discretized with the natural layer 
spacing d 1=3 12A. The total number of layers in the z- 
direction is denoted by L z . We employ open boundary 
conditions in this direction, which allows topological de- 
fects to enter the sample from the top and bottom. The 
in-plane coordinates R M satisfy periodic boundary condi- 
tions and lie in a parallelogram spanned by the unit vec- 
tors e a = e x and ib = (e x + V3e y )/2 with area given by 
the number of vortices N times $0 / B. The unit of length 
is taken to be the coherence length £(T), with £(0)=17A, 
so that -ff C 2(0)=120T; the average inter- vortex distance 
is a$. The natural unit of time is t = £, 2 i]/eo and the 
time t is discretized with timestep dto] we take 5=0.05. 
We mention two technical aspects of importance. Con- 
siderable gain is reached in computational speed by gen- 
erating and storing the pinning forces (the "pinscape" , 
see Fig. 2) on a regular array of spacing £. The pinning 



force at a specific vortex location is then constructed by 
interpolating between the four closest grid-points. Sec- 
ond, in order to equilibrate disordered systems, we em- 
ployed simulated annealing, cooling from ~ 2.5T to T 
in 10 % decrements. Most of our simulations were done 
for 7 x 7 vortices in 50 layers, using 100,000 sweeps for 
equilibration and 100,000 for measurement. With the 
exceptions below we found only limited fluctuations be- 
tween different disorder configurations, due to the good 
self- averaging of the ~ 2500 vortex elements. Therefore, 
typically we averaged over only a few disorder realiza- 
tions. Finally the disorder strength was chosen to yield 
reasonable low temperature values for the critical current 
in terms of the depairing current jo: j c /jo ~ 0.005. 

We now proceed to investigate the transitions between 
the three proposed phases, the dislocation free Bragg 
Glass (BG), the dislocated Vortex Glass (VG), and the 
Vortex Liquid (VL). 




FIG. 3. The vortex wandering B{z) in the z-direction for 
different magnetic fields at T/T c = 0.7, 7=0.1, and e=0.2 
shows the Bragg to Vortex Glass transition. 



The BG-VG transition. In the dislocation-free 
Bragg Glass vortices fluctuate in the self-consistent cage 
of their neighbours (m) . Increasing the field brings vor- 
tices closer together, and increases the relative strength of 
the disorder. When the single vortex wandering reaches 
the Lindemann limit, vortices entangle and dislocations 
enter the sample |l4| , p8f . Our measured vortex line wan- 
dering, B(z) —< [R M (z, t) — R M (0, t)] 2 > fa%, is displayed 
in Fig. 3 for a typical field sweep. The curves for the 
five lowest fields group together and belong to the Bragg 
glass, whereas the curves for the higher fields are in the 
vortex glass phase. Our low field data in the BG in Fig. 3 
grow slowly with z, and are consistent with a power law 
with exponent less than unity. The presently attainable 
system sizes are, however, too small to see the asymp- 
totic regime, where B(z) ~ ln{z) M|. Increasing the 
field across the transition enhances B(z) abruptly and 
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instead of bending down, the curves now bend upward 
until boundary effects become important around z ~ L z . 

The vortex wandering in the VG phase reaches the 
Lindemann number cl at the entanglement length in the 
z-direction Le- From B(Le) = c| < O.f, we deduce 
that just above the transition, the entanglement length 
is of the order of 15 lattice spacings in the z-direction for 
the parameters of Fig. 3. Thus, our system size in the 
z-direction (50) is large enough to see entanglement, as 
evidenced by the snapshot just inside the VG phase in 
Fig. 4b, which shows two entangled "ring-exchanges" , or 
screw-dislocation loops, involving 2 and 9 vortices. On 
the other hand, Fig. 4a, taken on the BG side of the 
transition shows a perfect disentangled lattice. 
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FIG. 4. Snapshots of a typical configuration projected in 
the xy-plane in the a) the Bragg glass phase at B=6T and b) 
the Vortex Glass phase at B=7.2T. Parameters as for Fig. 3. 

The average height of the six principle structure fac- 
tor peaks in S(Q) =< (1/JV) exp(£Q • [R M CM) - 
'R. u (z, t)]) > is also indicative of the BG-VG transition. 
This is shown in Fig. 5, together with the wandering 
across the sample B(L Z ). In the BG phase, the structure 
factor first increases with field, due to the growing inter- 
action strength. Upon crossing the BG-VG transition, 
the structure factor peak height drops at the same point 
where the wandering increases dramatically. Frozen-in 
short range order still persists in the VG phase, however, 
as is clear from Fig. 4b. 

For different disorder realizations, the jump in B(L Z ) 
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FIG. 5. The magnetic field dependence at T/T c =0.7 of 
the averaged amplitude of the first 6 primary structure factor 
peaks (boxes), normalized to the number of vortices iV=49, 
and the wandering across the sample B(L Z ) (diamonds). 



and the corresponding drop of the structure factor peaks, 
occur at different magnetic fields. As a result, a disor- 
der average of the quantities we measure shows only a 
smooth field dependence. Although the sharp jump in 
the wandering for a single disorder realization is consis- 
tent with a first order transition, we speculate that in 
the thermodynamic limit only a broad feature survives, 
which mimics a continuous transition or a cross-over. 

VG-VL transition. This transition is defined by a 
profound change in the dynamical behavior of the system: 
the vortex glass can be seen as a frozen liquid. This 
is most clearly seen in the behavior of the time overlap 
correlator C(t) = exp(- < [R M (z,t) - R^(z,0)] 2 > /l 2 ), 
as plotted for ?=0.2a$ in the inset of Fig. 6. For times 
t 3> to it decays exponentially with a "plastic" relaxation 
time t p i , which is plotted as a function of temperature in 
Fig. 6. Upon cooling from the VL to the VG, the time 
t p i rapidly increases, signalling that vortices get trapped 
in the almost static cages of the surrounding vortices. 
The plastic time does not appear to be infinite in the 
VG, however. It saturates for temperatures below the 
VG-VL line, at a value around t p i ~ 2000io- Whether 
this is a finite size effect, or hints to the absence of a 
genuine phase transition with universal exponents, will 
be the subject of a future finite size study |2{|. 

As a function of temperature and field, the six pri- 
mary structure factor peaks decrease, but in the vortex- 
glass phase with frozen short range order, they remain 
finite and visible in our finite system and simulation 
time. As the vortex glass melts into the liquid, the pri- 
mary peaks dissolve into a ringlike structure. Transport 
measurements can also identify the VG-VL transition. 
Fisher et al. jlj| suggested that the flux flow resistivity 
at high temperature gives way to non-ohmic exponential 
creep below the glass transition temperature, a prediction 
which was confirmed experimentally [pOj, see, however, 
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FIG. 6. The relaxation time t p ;/(10 3 to) as a function of 
temperature from an exponential fit to the correlator C(t) 
illustrates the VG-VL transition. Here B—4T and the dashed 
line is a guide to the eye. The inset shows the correlation 
function C(t) as a function of time for T/T c = 0.91, .., 0.95 
(step 0.005, full lines, from top to bottom) and the fits (dashed 
lines) . 

also Ref. [ pl[ . Our numerical results on the IV curves 
will be presented elsewhere p9lj . 

The above diagnostics allow us to construct the phase 
diagram. We identified the predicted dislocation free 
Bragg Glass, entangled Vortex Glass, and Vortex Liquid 
phases in the B—T plane, as summarized in Fig. 1. We 
find that the BG-VG phase boundary is rather sensitive 
to the disorder, whereas the VG-VL transition is not. In 
fact, even from disorder realization to realization, the lo- 
cation of the BG-VG boundary shows detectable fluctua- 
tions. This sensitivity is a well known experimental fact: 
Safar et al. report transition fields from 3T to 10T 
in samples with minute macroscopic differences, whereas 



the data of Ref. |17| are at even lower fields. Experimen- 
tally, the BG-VG transition appears as a sharp jump in 
the magnetization only in BSSCO, while in YBCO and 
Thallium based materials the transition is rather smooth. 
The reported numerical fluctuations seem to be related 
to this smearing of the BG-VG transition in the latter 
compounds. In Fig. 1 we also display the experimental 
data of Deligiannis et al. []l9| and Safar et al. (2^|. We 
chose a disorder value to make the transition field fall 
between the two sets of data. The correspondence be- 
tween the simulations and experimental data is remark- 
able, especially for the VG-VL transition. At the highest 
temperatures T/T c ~ 0.95, the BG-VG phase boundary 
bends downward. The resolution of our data is insuffi- 
cient at the moment to decide, whether the BG-VG line 
ends somewhere on the melting line, as the BSSCO ||] 
and some YBCO data VM indicates, or whether there is 



a sliver of the VG phase intruding all the way down to 
zero fields, preempting a direct BG-VL transition, as the 
Thallium Jl§ and other YBCO data 0Jl9) suggests. 
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